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In paper III of this series we proposed a scenario of superconductivity driven by hole "undressing" 
that involved a complete redistribution of the occupation of single particle energy levels: the holes 
near the top of the band were proposed to all condense to the bottom of the band. Here we consider a 
less drastic redistribution involving electrons with a definite spin chirality and show that it is in fact 
energetically favored by the Coulomb exchange matrix element J over the scenario proposed earlier. 
It is shown that spin splitting with chiral states reduces the Coulomb repulsion and hence that the 
Coulomb repulsion promotes spin splitting. Superconductors are proposed to possess a spin-split 
hole 'core' at the bottom of the electronic conduction band in addition to a spin-split Fermi surface. 
The new scenario leads naturally to the existence of a spin current in the superconducting state and 
is consistent with the Spin Meissner effect and negative charge expulsion discussed earlier within 
the theory of hole superconductivity. 

PACS numbers: 



I. INTRODUCTION 

The theory of hole superconductivity proposes that su- 
perconductivity can only occur when the Fermi level is 
near the top of an electronic energy band and that it is 
driven by "undressing" of charge carriers that are hole- 
like in the normal state and become electron-like in the 
superconducting state[l|. The microscopic interaction 
proposed to drive this physics is a correlated hopping 
term At 2} arising from an off-diagonal element of the 
Coulomb interaction in a tight binding representation Q 
in the presence of orbital relaxation [4J, that is attractive 
for states near the top of a band and particularly strong 
when the ions are negatively charged] 5], and that makes 
the band increasingly more narrow [6| and incoherent [7j 
as the electronic occupation increases. 

Analysis of the consequences of this theory for the 
electrodynamics of superconductors led to the conclu- 
sion that drastic modification of the conventional Lon- 
don electrodynamics is required: a macroscopically in- 
homogeneous charge distribution is proposed to exist in 
superconductors, with an internal electric field point- 
ing towards the surface[8|- Associated with it is a 
spin current flowing within a London penetration depth 
of the surfaced. The process by which the charge 
inhomo geneit v |1G| and the spin current [ll| are generated 
in the transition from the normal to the superconducting 
state provides a 'dynamical' explanation of the Meissner 
effect [12[, which does not exist within the conventional 
London-BCS theory [3. 

We have also shown in previous work that the inter- 
action At yields a lower energy in the presence of a spin 
current in the superconducting state jl5|. However, the 
microscopic origin of the spin current was not clarified 
in that work. Furthermore, in work unrelated to super- 
conductivity we have shown[l|| that another off-diagonal 
matrix element of the Coulomb interaction, J, favors a 



new type of Fermi surface instability of metals, a 'spin- 
split state' where Fermi surfaces of opposite spin elec- 
trons shift relative to each other and a spin current de- 
velops, with that state being particularly favored when a 
band is nearly half-filled. In related work, Wu and Zhang 
have recently discussed dynamic generation of spin orbit 
coupling within a Fermi liquid framework [j~7j and cast the 
problem within the general framework of Pomeranchuk- 
type instabilities in the spin channel[18]. 



The purpose of this paper is to propose a new sce- 
nario for the generation of a spin current within the the- 
ory of hole superconductivity. In the previous paper in 
this series we proposed that holes 'undress' from the 
electron-ion interaction by migrating from the top to the 
bottom of the band in the transition to superconductiv- 
ity, and argued that this provides an explanation for the 
negative charge expulsion and resulting internal electric 
field predicted by the electrodynamics equations. How- 
ever, that scenario did not address the question of how 
the spin current develops. Here we show that a modifica- 
tion of that scenario does. Namely, it is not holes of both 
spin orientations that migrate to the bottom of the band, 
but only holes of one definite chirality. The chirality di- 
rection is determined by the Dirac spin-orbit coupling 
in the presence of the electric field generated by charge 
expulsion, which acts in an analogous way as a 'seed' 
magnetic field would act in a ferromagnet to determine 
the preferred broken-symmetry state. A driving force for 
this effect is shown to be again the off-diagonal matrix 
element J, which lowers the energy for a superconduc- 
tor in this scenario when the band is almost full rather 
than near half-filling as was the case for the spin-split 
'normal' state proposed in Ref . (l6j . Furthermore, a key 
role is proposed to be played by the ordinary Coulomb 
interaction and in particular its long-range nature. 
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FIG. 1: The diatomic molecule as a microcosm for an energy 
band. When the electron-electron interaction dominates, the 
occupation of single-particle energy levels changes, as shown 
on the right side of the figure: the middle diagram represents 
a ferromagnet and the upper diagram a superconductor. 
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FIG. 2: The plane z = is the boundary of the supercon- 
ductor. An electric field in the interior of the superconductor 
points in the +z direction. There is excess negative charge 
within distance Al from the surface of the superconductor 
that was expelled from its interior, and there is a spin current 
near the surface. The long and short perpendicular arrows 
denote the direction of the momentum and the carrier's spin. 



II. THE NEW SCENARIO 
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We start by reproducing a figure from the earlier pa- 
per in this series 19] in Fig. 1, with its caption, depicting 
the states of a diatomic molecule with 3 electrons. We 
argued in Ref.[lj| that the Coulomb matrix element At 
lowers the energy when two electrons are in an antibond- 
ing orbital and hence will promote electrons from lower 
to higher orbitals. Then we extrapolated to a lattice sys- 
tem and argued that all the electrons near the bottom of 
the band will be promoted to the top of the band, as was 
depicted in Fig. 6 of ref. [19| ■ This would indeed be the 
case if the effective hopping amplitude as function of the 
number of electrons per site n e 



t(n e ) = to — n e At 



(1) 



were to change sign for n e approaching 2, the full band. 
That this is not a totally implausible assumption is 
shown in Appendix A. However, we believe neverthe- 
less that it is too stringent a requirement. In fact, for 
the high T c cuprates comparison with experiments led us 
to conclude[20] that while t(n e ) becomes very small as 
n e — ^ 2 it does not change sign. 

Nevertheless, for the diatomic molecule the right panel 
of Fig. 1 also describes schematically the scenario pro- 
posed in this paper. Which is that half of the electrons 
in states near the bottom of the band get promoted to 
near the top of the band when a system becomes super- 
conducting. 



III. SPIN CURRENT AND DIRAC PHYSICS 

Within the theory of hole superconductivity the su- 
perconductor expels negative charge from the interior to- 
wards the surface [l^, resulting in an excess of negative 
charge density within a London penetration depth Al of 



' Am e c\L 



(2) 



with n s the superfluid carrier density and Al the London 
penetration depth given by 



Al = ( 



m e c 
47m, e 2 



1/2 



(3) 



In the same region near the surface there is a spin current 
flowing parallel to the surface, with carrier speed given 



bylUl 



4m P A; 



(4) 



The magnitude of v® can be understood from the fact 
that it gives rise to an gul ar momentum h/2 for carriers 
in orbits of radius 2Al|ii|, [Hj], and orbits of radius 2Al 
are required to understand the Meissner effect [TTI. Il3j . 
The carrier spin a points parallel to the surface and per- 
pendicular to its velocity v according to the relation^ 

a = v x n (5) 

with n the normal to the surface pointing outwards. 

Figure 2 shows schematically a superconductor in the 
region z < bounded by the z — plane. The Dirac 
spin-orbit interaction of carriers near the surface in an 
electric field E pointing towards the surface is[2lj 



Hs.n. — 



: <7 • (E x p) 



(6) 



where E is the electric field and p is the momentum. 
The magnitude of the electric field in the superconduc- 
tor near the surface according to the theory of hole 
superconductivity Q is E m given by 



he 



(7) 
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Hence Eq. (6) becomes, with p = fik 



h 2 

2m, 



kqo( 



2m e c 



qo)a ■ (k x z) 



where we have defined 



1 

2A? 



(8) 



(9) 



A spin-orbit term of the form Eq. (8) is commonly re- 
ferred to as a Rashba term_22j and the corresponding 
spin-split bands as Rashba bands. 
The quantity 



r„ = 



2m P c 



(10) 



(proportional to the Compton wavelength) may be called 
the 'quantum electron radius' (as opposed to the 'clas- 
sical electron radius' r c = e 2 /m e c 2 ) : a mass m e or- 
biting at speed c with radius r q has angular momen- 
tum h/2, the electron spin, and furthermore the quan- 
tum confinement energy of a mass m e in a distance r q is 
h 2 /(2m e r 2 ) = 2m e c 2 . For a London penetration depth 
Xl = 160A (Al) the quantity in brackets in Eq. (8) is 



r q qo 



6.0 x 10 



-6 



(11) 



that sets the scale of the spin-orbit symmetry-breaking 
field arising from Dirac single-particle physics. The spin 
current speed resulting from Eq. (8) is 



Dirac 



2m f 



-qo(r q q ) = (r q q )v° 



(12) 



and the favored spin direction is given by Eq. (5) . 

Thus the situation is similar to that of of fcrromag- 
netism. The magnetic dipole interaction between elec- 
tronic magnetic moments of magnitude fiB (/is=Bohr 
magneton) can provide a qualitative explanation of fer- 
romagnetism (or antiferromagnetism) but its strength is 
several orders of magnitude too small to explain the ob- 
served ferromagnetism in solids. Hciscnbcrg correctly 
surmised that the true explanation of ferromagnetism 
lies in the much stronger Coulomb interaction between 
electrons. Similarly here, the spin-orbit interaction aris- 
ing from the Dirac equation gives rise to a spin current 
speed which is 10 -5 times smaller than what the electro- 
dynamic equations of the superconductor tell us it should 
be[9j, namely Eq. (4). This is because the ultimate ori- 
gin of the spin current is again the Coulomb interaction, 
just as in the case of ferromagnetism. 



IV. ORIGIN OF THE SPIN CURRENT AND 
ELECTRIC FIELD 

In our earlier work on metallic ferromagnetism |23j as 
well as the spin-split state[l6j], we argued that these col- 
lective effects are driven by a Coulomb matrix element J 
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FIG. 3: Schematic image of the Brillouin zone occupation 
in the superconducting state for an s-like electronic energy 
band (bottom of the band at the F-point). Depicted is a 45° 
plane that goes through the points (0, 0, 0) and (n, n, it) in 
the Brillouin zone projected onto the (k x ,k y ) plane. The k z 
axis points out of the paper. A circle of radius qo = 1/2Al 
in the (k x ,k y ) plane centered at (0,0) is singly occupied by 
electrons of chirality a ■ (k x z) — — 1. The arrows near 
the center indicate the direction of (k x ,k y ) (radial arrow) 
and of the spin of the corresponding electron (attached tan- 
gential arrow). There is an excess of electrons of chirality 
3 • (k x z) = 1 at the Fermi surface in the antibonding regions 
close to (k x , k y ,k z ) — (in, ±7T, ±7t), (their spin orientation is 
indicated by the arrows). The resulting regions with singly- 
occupied electrons are denoted by the hatched areas, with 
45° hatching (-45° hatching) in the antibonding (bonding) 
regions respectively. The region hatched by the square pat- 
tern is doubly occupied with electrons, and the white regions 
near (k x ,k y ,k z ) = (±7r, ±7t, ±7t) are doubly occupied with 
holes. 



( "bond-charge repulsion" [2J, [25[ ) that favors the Fermi 
surfaces of opposite spin electrons to occupy regions of 
the Brillouin zone of opposite bonding character (i.e. one 
bonding, one antibonding). We invoke the same physics 
here. 

Consider for definiteness a three-dimensional cubic lat- 
tice with an s— like band close to full. In real space, the 
electric field originating in the charge expulsion points 
along the +z axis as shown in Fig. 2. According to 
the Dirac interaction Eq. (8), electrons with spin ori- 
entation satisfying a ■ (k X z) = +1 arc slightly fa- 
vored over those with opposite chirality. We propose 
that in order to lower the bond-charge Coulomb repul- 
sion energy, half of the electrons in a small region near 
k ~ (bonding character) will migrate to the region 
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(k x , k y ,k z ) ~ (±7r, ±7r, ±7r). Which half of the electrons? 
Those with chirality, a ■ (k x z) = 1. The resulting situ- 
ation is shown schematically in Figure 3. 

As in Ref.pjjj, we argue that the electrons expelled 
from the center of the Brillouin zone occupying the more 
costly electron-ion energy states near the corners of the 
Brillouin zone represent the expelled electrons from the 
interior of the superconductor that reside near the real 
space surface. Their chirality corresponds to the direc- 
tion depicted in Fig. 2 predicted by the Spin Meissner 
effect |nj]. 

As Fig. 3 shows, we end up with a new 'Fermi surface' 
deep in the interior of the Brillouin zone, of radius qo = 
1/(2Al). The carriers in that region of single spin chiral- 
ity occupy long- wavelength states (q^ 1 >> lattice spac- 
ing) that don't 'see' the discrete nature of the electron-ion 
potential, hence they are 'undressed' from the electron- 
ion interaction [2^| and behave as free electrons (or holes). 
They 'drive' the behaviour of the entire superfluid con- 
densate in the superconducting state. The process by 
which the holes near {k x ,k y ,k z ) = (±7r, ±7r, ±7r) mi- 
grate to k ~ [26[ in the transition from the normal to 
the superconducting state and acquire chirality is shown 
schematically in Fig. 4. It corresponds in real space to 
the orbit expansion from an orbit of atomic dimension 
(corresponding to a nearly filled band) to an orbit of ra- 
dius 2A_l, as predicted by the Spin Meissner effect [11 1 
and shown schematically in Figure 5. Note that all elec- 
trons in the Fermi sea are affected by this process, which 
is consistent with the physics of the Meissner effect p7l|. 
In contrast, within BCS theory only a small fraction of 
electrons near the Fermi surface are affected in the transi- 
tion to superconductivity. The orbit expansion is neatly 
represented by the change in the Larmor diamagnetic 
susceptibility 



XLa 



n s e 2 
-(a) = — - -a 



Am e c 2 



(13) 



for orbits of radius a from 



k tp 



corresponding to the 




FIG. 4: The figure shows a pair of holes (k f, — k J.) with 
spin orientation k x z dropping from the Fermi level to the 
bottom of the electronic energy band. In the process k and 
k rotate around the k z direction and Ikl shrinks. 
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FIG. 5: The figure shows a superconducting cylinder viewed 
from the top. Electronic orbits in the normal state have ra- 
dius kp 1 , of order of the ionic lattice spacing, and electronic 
orbits don't overlap. In the transition to superconductivity 
the orbits expand to radius 2Al, several hundreds Angstrom, 
and they become highly overlapping. Only orbits with nor- 
mal parallel to the cylinder axis are shown. Electrons with 
spin pointing out of (into) the paper traverse the orbits in 
clockwise (counterclockwise) direction in the superconductor. 



1 2 , s 
XLandau — -^fJ-B9\ e F) — XLa 



ria^kp 1 ) (14) 



(with g{tF) = 3n s /2eF the density of states) to the 
London susceptibility in the superconducting state cor- 
responding to Larmor orbits of radius a — 2Al = Qq^IHI 



XLo 



1 

47T 



XLa 



.(a = q 1 ). 



(15) 



The large orbits depicted on the right part of Fig. 5 are 
highly overlapping, and for that reason phase coherence 
is a must to avoid collisions between electrons that would 
be very costly in Coulomb energy [2^. This is why a state 
with the occupations depicted in Fig. 3 cannot exist as a 
normal metallic state, but requires the long range phase 
coherence characteristic of a BCS-like wavefunction. 



V. COULOMB MATRIX ELEMENTS 

Figure 6 shows qualitatively the form of the wavefunc- 
tion for electrons near the bottom and near the top of a 
band. In the region between the ions the charge density 
is large for electrons near the bottom of the band (bond- 
ing states) and small for electrons near the top of the 
band (antibonding states). This is of course the reason 
for the name: the large charge density between the ions 
of the bonding electrons give rise to an effective attrac- 
tion between the ions, thus bonding the crystal together. 
When a band becomes almost full with many electrons in 
antibonding states the crystal becomes unstable because 
of the lack of interstitial charge density in the occupied 
antibonding states. 
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FIG. 6: Qualitative nature of electron energy states near the 
bottom and near the top of an electronic energy band. For the 
states near the bottom (bonding states) the charge density is 
large in the region between the ions, for the states near the 
top (antibonding states) it is small 



Figure 6 suggests that the Coulomb repulsion energy 
will be lowered if some electrons occupy bonding states 
and some occupy antibonding states, compared to the 
case where all electrons occupy bonding state, because 
the bond charge revulsion^2/&, |25|| between electrons will 
be reduced. We have argued that this effe ct p lays a fun- 
damental role in metallic ferromagnetism[]23j and in the 
spin-split state[l6j], and we propose here that it also plays 
a fundamental role in superconductivity, by giving rise to 
the Brillouin zone occupation depicted in Fig. 3. 

Consider the nearest neighbor matrix elements of the 
Coulomb interaction between s-wave like Wannier states 
<Pi 

(ij\kl) = J d\d\'rt{r)^(r')^-^w{r') Vk {r) (16) 

giving rise to U = V = J = 

J' = At = all positive. The single particle 

hopping and the 'hybrid' matrix element At give rise to 
the kinetic part of the HamiltonianQ 

H k in = - y^Jto - At(7tj,_ <r +n j - (T )][cl (T Cj <T +h.c.]. (17) 

The nearest neighbor exchange interaction J gives for 
antiparallel spins 

Hj = J J2 c lcjA-^- ( 18 ) 

and the pair hopping J' 

Within a mean field decoupling of these terms (23l. |29j . 
the effective hopping amplitude for an electron of spin a 
is 

t = t - n e At - ( J + J')I- a (20) 

with 

n e =< nit > + < nn > ( 21a ) 



the site occupation and 

I-a =< 4,-<rCj,-° > ( 21b ) 

the bond occupation for the opposite spin orientation. 

The matrix elements J and J' involve overlaps of two 
Wannier orbitals and hence are expected to be negligible 
beyond nearest neighbors. The single particle hopping 

as well as the hybrid matrix element At involve a sin- 
gle orbital overlap and thus could be more long-ranged. 
Hence we assume there is also a next-nearest-neighbor 
hopping amplitude 

t 2 = t 02 - n e At 2 . (22) 

As shown in Appendix A we expect that the ratio At 2 /to2 
should be larger than At/t , and we will assume that the 
parameters satisfy to — n e At > and to2— n e Ati < 0. For 
a two-dimensional square lattice the dispersion is then 

e^k = -2t nn (cos(k x ) + cos{k v )) 

+2t nnn cos(k x )cos(k y ) (23a) 

t nn = t - n e At — ( J + J')I- a (23b) 

tnnn = n e At 2 - t 02 (23c) 

with t nn > 0, t nnn > 0. The bond charge occupation Jo- 
depends on temperature and will increase as the temper- 
ature decreases. We have discussed in earlier work how 
transitions to ferromagnetic [23| and spin-split 16;] states 
can occur when I a reaches a critical value for a half-filled 
band. Similarly we envision here that for a band that is 
almost full as the temperature is lowered t nn decreases 
relative to t nnn due to the increase in the bond charge 
occupation I a and consequent increase in the electron- 
electron bond charge repulsion represented by J and J' 
and a point will be reached where it becomes favorable to 
transfer electrons from k ~ to k ~ ±7r. This is shown 
schematically in Figure 7. 

Below the critical temperature, we expect that the ex- 
pectation value of the bond charge will break parity, as 
in the spin-split state[l^]. For a given k direction, the 
bond charge corresponding to a ■ (k X z) = — 1 decreases 
and the opposite one increases. It can be seen from Eq. 
(23b) that when 7_ CT decreases the nearest neighbor hop- 
ping amplitude increases and thus the band energy corre- 
sponding to the opposite spin orientation increases. The 
spin-split bands will be given by 

h 2 k 2 h 2 

(•ka = kq a ■ (kx z) (24) 

2m e 2m e 

and holes of chirality a ■ (k X z) = 1 will condense to 
the bottom of the lower spin-split band in a circle of 
radius q . A superconducting gap of magnitude 40] Aj, = 
h 2 kqQ/2m e opens up at the chemical potential, as shown 
schematically in Figure 8. 
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FIG. 7: Evolution of the single particle dispersion (Eq. (23a) 
versus temperature for a two-dimensional square lattice along 
the line k y — k x . The horizontal lines denote the position of 
the chemical potential for band occupation 0.1 holes per site. 
The parameters for the top, medium and bottom panels are 
tnn = 1.5, 0.3 and 0.16 with t nnn = 0.5. The bottom panel 
shows the threshold where holes at the center of the Brillouin 
zone start to become energetically favored. 



VI. COULOMB INTERACTION ANALYSIS 

The Coulomb interactions resulting from the Coulomb 
matrix elements Eq. (16) are given in momentum space 

by 

H uv = ^ Yj T/ ( c l)4+q<T C L-q ( T' c k'<T'C k(T (25a) 

kk'qcrcr' 



Hj= 2N ^ J(k ~ k ' + q ) c k+ qff c k'- qCT ' c kV'C kCT 

kk'qerer' 

(25b) 

Hj ' = ^N £ J/ ( k + k ')Ck +qCT 4'- qCT 'Ckv<c kCT (25c) 

kk'qo'O'' 



-2 2 

k x /rr 

FIG. 8: Spin-split bands below T c given by Eq. (24). fcFisthc 
Fermi wavevector, and /i = h 2 k F /2m e . The dot-dashed line 
denotes the energy dispersion relation in the normal state. 
The circle of radius go in the lower band, corresponding to 
a ■ (k x z) = 1, is unoccupied by electrons, i.e. is occupied 
by holes. Note that in the region kF — go < k < kF only 
electrons of chirality a ■ (k x z) = 1 exist. 



H " = 2N S [Ai(k)+At(k')+Mk + q) + 

kk'qcrcr' 

At(k' - q)] x 4 +q a c L- qCT ' c kvc kCT (25d) 
with the interaction amplitudes given by 

V(q) = U + V^e 1 * 5 (26a) 

<5 



J(q) = jJ2e^ s (26b) 

<5 



J'(q) = J'5>"^ (26c) 

s 



At(q) = At^e 1 * 8 (26d) 

5 

with (5 connecting a site to its nearest neighbors. These 
interactions are maximally positive for q = and maxi- 
mally negative for q at the corners of the Brillouin zone 
q = (±7r, ±7r, ±7r) = ±7?. Assuming that superconduc- 
tivity is driven by a reduction of the Coulomb interaction 
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between electrons, let us consider the conditions that lead 
to such reduction for the various terms. 

The self-energy of the electrons is obtained from the 
expressions Eq. (25) for q = 0. From Eq. (25d) it can be 
seen that At gives negative contributions when k and k' 
are close to the corners of the Brillouin zone[3(| IHj], thus 
favoring bands that are almost full and have hole conduc- 
tion in the normal state, in accordance with Chapnik's 
rule(32|. J(k — k') will give negative contributions when 
k — k' ~ ±7? , so it favors simultaneous occupation of the 
center (k ~ 0) and the corners of the Brillouin zone, as in 
the spin-split scenario discussed in the previous section. 
Similarly J'(k + k') will give negative contributions in 
such a case. V(q = 0) is uniformly repulsive indepen- 
dent of the occupation of the Brillouin zone. 

Next let us considering the scattering terms in the BCS 
channel: 

1 



tjBCS 
H UV 



N 



E 

kk' 



v(k - k')4 



kf -ki 



C-k'4.Ck't 



(27a) 



H J CS = jf E J ( k + *')<iAn c -*w (27b) 

kk' 

H f'° S = Jj E ./'(0)<W_ ka c k , t (27c) 



kk' 



Hl? S = i^[At(k) + At(-k) + At(k') + At(-k')] 



kk' 

"kt C -ki C 



-k'j.c k 't (27d) 

The At terms give negative contribution for scattering 
when k and k' are near the corners of the Brillouin zone 
and thus give rise to hole pairing, as previously discussed 
extensively [3. l20j. J gives negative contributions when 
k + k' is close to the corners of the Brillouin zone, hence 
one pair in the 'inner' Fermi surface and the other pair in 
the 'outer' Fermi surface in the spin-split state depicted 
in Fig. 3. J' gives repulsive contribution regardless of 
the locations of k and k'. 

Most significantly, V(k — k') will be least repulsive for 
k k' close to ±7?, since the nearest neighbor interac- 
tion contribution becomes attractive (Eq. (26a)). This 
fact has been emphasized by Mattis in the context of 
an electrostatic mechanism for superconductivity [30j . In 
the present scenario, it will lower the Coulomb energy 
for pair scattering between the inner and outer Fermi 
surfaces depicted in Fig. (3). 

In summary, we conclude from this analysis that the 
proposed scenario, with an almost full band in the nor- 
mal state and spin-splitting as depicted in Fig. 3 in the 
superconducting state, appears to be optimal to reduce 
the various Coulomb interaction terms arising from the 
matrix elements of the Coulomb interaction in the Bloch 
states of the electrons. As discussed elsewhere [33T - I351 ] . 
the off-diagonal matrix elements of the Coulomb inter- 
action are augmented by orbital relaxation effects which 
are largest for negatively charged ions. 



VII. BCS WAVEFUNCTION AND COULOMB 
INTERACTION 

We consider a BCS-likc wavefunction of the form 



l*>=Il( Uk 



U kC kt cL k4 



)(«-!< 



«-kcVk) (28) 



where the product is over half the Brillouin zone (e.g. the 
subspace k x > 0) and where we denote by (k f , — k J.) 
pairs of favorable chirality (<?• (k x z) = 1) and by (— k j" 
, k J,) pairs of unfavorable chirality. Quite generally, spin- 
splitting will be described by the condition 



M 2 ^ Kkl 



(29) 



in some region(s) of the Brillouin zone. As discussed ear- 
lier, this parity breaking would be caused just by the 
Dirac spin-orbit interaction in the presence of the elec- 
tric field generated by charge expulsion, but the effect 
will be greatly augmented by the Coulomb interaction, 
increasing the characteristic wavevector from go x ( r qQo) 
to g . 

Consider scattering processes of the form (k f, k I 
) —> (— k f, k l) coupling the two Rashba bands. The 
associated expectation value for the Coulomb matrix el- 
ement V(k - k') = V(2k) is 



V(k - k')M k u k M!l k u_ k 



(30) 



so the interaction is repulsive both for states near the 
original Fermi surface (k ~ ±7?) and near the new small 
Fermi surface (fc ~ go ~ 0). Spin splitting will reduce 
this repulsive energy because UkW k and uL k v_ k will not 
be simultaneously large for the same k. The same is 
true for other scattering processes coupling both Rashba 
bands. Thus, in the superconducting state V(k) provides 
a strong driving force for separation of the Fermi surfaces 
of opposite chiralities initially split by a tiny amount due 
to Dirac spin-orbit coupling. 

Furthermore for scattering within the same Rashba 
band V(k — k') can give an attractive interaction and 
corresponding energy lowering. In computing the ex- 
pectation value of the Coulomb interaction we need to 
include the overlap matrix elements M(k, k') between 
Rashba spinors, since the orientation of the spin changes 
with wavevector. For the small Fermi surface the Rashba 
spinors are 



1 ( 1 

V2 \ Vl+fcf 



(31) 



and the overlap matrix elements for scattering from an 
initial state k = (go,0) to k' = (k x ,±yjq^ — k%) add up 
to 



Mi + M 2 = 



k x (k x + g ) 



?o 2 



(32) 
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so that the interaction is attractive when the angle be- 
tween k and k' is between 90° and 180° (—qo < k x < 0). 
Similarly, in scattering processes between the hole pock- 
ets near ±7? involving the same Rashba band, the relative 
angle between spinors will be (approximately) either 90° 
or 180° (see Fig. 3), hence will give vanishing overlap 
matrix elements (corresponding to taking k x = and 
k x = —qo in Eq. (32)) and suppress the repulsion due to 
y(k-k'). 

For the interaction At similar effects occur, so in par- 
ticular the attractive interaction for scattering between 
hole pockets near ±7? will be weakened by the spin split- 
ting. However, it should be remembered that the magni- 
tude of At is much smaller than that arising from diago- 
nal elements of the Coulomb interaction, e.g. the nearest 
neighbor repulsion V , so the net effect certainly should 
be a tendency to enhance spin splitting. For scattering 
between states near the small Fermi surface the inter- 
action At is strongly repulsive so the beneficial effect of 
spin splitting also applies to it. 

For scattering processes that do not change (k x ,k y ) 
however this supression of the Coulomb repulsion due to 
spin splitting does not take place. Consider scattering 
processes involving the same k— components in the xy 
plane and opposite k z values 



k — {k x , k y , k z ) 



{kx i ky , k z ) 



(33a) 



(33b) 



The expectation value of the Coulomb interaction involv- 
ing scattering processes (k f, — k \) — > (k f, — k \) (in the 
same Rashba band) is 

V{k - k> k u k ^ k = V{2k z )u k vtulv i[ (34) 

and will be large and repulsive for small k z since both 
w k u k and u ^ v k wm be large for (k x ,k y ) on the 'small' 
Fermi surface. To reduce the Coulomb repulsion in 
this case we suggest that the Bloch wave fermion op- 
erators in Eq. (28) should be replaced by the following 
operators [1^: 



V2 



J 

--k| 



-k| 



V2 



(35a) 



(35b) 



and similarly for the opposite chirality operators. The 
expectation value corresponding to Eq. (34) with the 
wavefunction 

I* >= I]> k + w k2 k T 5 -ki)(w-k + w- k 5l kt 4i) ( 36 ) 

k 

is now zero (except for the particular case k, = 0), indi- 
cating that this basis will lower the energy [3g. In ad- 
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FIG. 9: Wave function amplitudes and charge densities along 
the z direction for the standing waves associated with the op- 
erators (full lines) and cl kCT (dashed lines) , assuming the 
Bloch states are plane waves. Note that the charge densities 
avoid each other reducing the Coulomb repulsion. 



dition to reducing the Coulomb repulsion for scatter- 
ing processes as seen above, the Coulomb self-energy 
for pairs (c k ^cl k ^) as well as for pairs (cL k *c k |) is re- 
duced compared to ordinary Cooper pairs, from V r (0)/f2 
to (7(0) - V(2k z )/2)/n) (fl ^volume). This is because 
these operators create standing waves in the z direction 
that are 90° out of phase between k and — k, as shown 
schematically in Fig. 9. 

The z axis direction corresponds to the direction of 
the electric field, which is normal to the surface of the 
superconductor. Thus it is reasonable to use in the BCS 
wavefunction Eq. (36) standing waves in the z direction 
because a superconductor in the ground state or in the 
Meissner state will carry charge and/or spin current in 
direction parallel to the surface but not in direction per- 
pendicular to the surface. For bodies of more general 
shape (e.g. spherical or cylindrical) the direction of the 
normal to the surface will change with position and the 
coordinate system needs to be changed accordingly. 

Note that the physics discussed here is in some sense 
opposite to BCS. BCS singled out (k f, k !) pairs be- 
cause they provide the largest phase space for scatter- 
ing onto other (k' |, k' \) pairs on the Fermi surface, 
which is good provided the scattering potential is attrac- 
tive (electron- phonon interaction) . Here we have instead 
a repulsive interaction and have argued that spin split- 
ting will reduce such scattering processes on the Fermi 
surface for large momentum transfers (eg (k f, k I) —> 
(— k f, k l) and even give negative contribution in a range 
due to the spin overlap matrix elements. 

We point out that it has recently been shown by Cap- 
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pelluti, Grimaldi and Marsiglio[37j that the tendency 
to superconductivity in systems with Rashba spin-orbit 
interaction will be greatly enhanced when the Fermi 
wavevector is comparable to the Rashba wavevector be- 
cause of an enhanced density of states. That is precisely 
the situation here (Fig. 8), suggesting that the scattering 
processes discussed involving the 'small' Fermi surface 
will have a dominant effect. 

What is the extent of the inner hole pocket in the k z 
direction? One can imagine various possibilities for the 
shape of the small Fermi surface, e.g. (i) a sphere of 
radius go > (ii) a prolate or oblate spheroid with axis along 
the k z direction, (iii) a cylinder of radius go extending 
all the way to the edges of the Brillouin zone, or (iv) 
a horizontal toroid of inner radius and outer radius 
go 37]. The answer is not clear at this point. It is also 
possible that it depends on the position of the point in 
real space relative to the surface, which would require a 
Bogoliubov-De-Gennes description. 

The following argument suggests that the inner hole 
pocket could be a cylinder extending to k z — ±w/a. The 
cost of expelling electrons from the region (k x ,k y ) ~ 
(0, 0) to the corners of the Brillouin zone gets progres- 
sively smaller as k z increases as far as the electron-ion 
and kinetic energy costs are concerned. The Coulomb 
energy cost for V(k — k') for processes between the inner 
and outer Fermi surfaces progressively increases. How- 
ever that cost could be offset if the gap function Ak 
were to change sign between k ~ (jr, tt, tt) and k' = 
(0,0, k' z ), rendering a repulsive contribution attractive. 
In fact, in the model of hole superconductivity originally 
considered [1, [2(| it was found that Afe changes sign for 
£k going down into the band where Ai(k) becomes re- 
pulsive, and this provides additional enhancement to T c 
('spatial pseudopotential effect '0, HH). In the current 
scenario we would expect Ak' to change sign for large 
k z but not for k z — > since in that case V(k — k') is 
attractive. 

Finally, we point out a peculiar feature in our earlier 
work on the model of hole superconductivity @, H3|- In a 
hole representation we had, as in usual BCS 



M 2 = 2(i- 

for the hole occupation, with 



and the gap function depended linearly on energy 



A k = A m (- 



D/2 



(37) 



(38) 



(39) 



with D the bandwidth and A m , c constants. For large 
ek, meaning states deep inside the Fermi sea 



(40) 
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FIG. 10: BCS coherence factors for the correlated hopping 
Hamiltonian of Ref . [39l] . in hole representation. The param- 
eters used correspond to Figure 1(a) of Ref. 39], with hole 
occupation nh = 0.17 and T c = 18K. D is the bandwidth, 
and the chemical potential /i is indicated by the dotted line. 
Full and dashed lines are Vk and ut in hole representation, the 
hole occupation is 2|«k| 2 . Note that the amplitude for hole 
occupation Vk is non-zero all the way down to the bottom of 
the electronic band (e.k/{D/2) = 0.5). 



and 



1, 1 

2 



D/2> 



=)>0 



(41) 



indicating that the hole occupation does not go to zero 
even for states deep inside the Fermi sea, unlike conven- 
tional BCS theory. An example is shown in Fig. 10, and 
this behavior is generic for this model. This could have 
been seen as a precursor of the scenario discussed in this 
paper. 



VIII. THE KEY TO SUPERCONDUCTIVITY 

The arguments in the previous section indicate that: if 
there is a 'small' Fermi surface deep inside an almost full 
band, it better be spin-split, i.e. corresponding to spins 
of one chirality only. But the arguments don't prove that 
there is a small Fermi surface. The arguments involving 
J are suggestive but perhaps unconvincing for the case 
of nearly-free-electron metals that become superconduct- 
ing. Furthermore, the arguments given so far do not give 
indication for why the degree of spin splitting in the su- 
perconductor should be go. 

Consider electron states deep in an almost full band, 
with wavevectors of order go. The Coulomb interaction 
certainly has a matrix element connecting the state (k f 
, -k I) with the state (-k f, k I): 

1 47re 2 

< -k t, k | \V c \k t, -k |>= V c (2k) = - _ (42) 

We have not included in Eq. (42) the Thomas Fermi 
wavevector to screen the long-range Coulomb interaction 
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as is usually done. We argue that this interaction be- 
tween states deep in the Fermi sea cannot be Thomas- 
Fermi screened because the electronic states needed to 
be reorganized in order to screen the interaction are oc- 
cupied by other electrons and thus not available. Thus 
these matrix elements become very large when k becomes 
very small. 

The reader may argue that because all those states 
are occupied in the normal state Fermi sea, the Coulomb 
interaction is unable to scatter (k f , k 1) into (— k "J" 
, k 1) because of the Pauli exclusion principle, hence that 
such processes are irrelevant. 

However, there are other states with such wavevec- 
tors and spin orientations: the negative energy states 
in Dirac's theory of the electron. 

If the pair (k "f, k 1) makes a transition to the top of 
the band, it opens up room for a pair of electrons in the 
Dirac sea of negative energy states to make a transition 
to that state. 

To create a Dirac electron-hole pair costs an enormous 
energy, 2m e c 2 . However, the Coulomb matrix element 
mixing a pair (k f, k |) with another pair (— k f,k I) 
and |k| = qo has magnitude 



V(2q ) 



n 



(43) 



with N s /2 the number of superfluid pairs. 

Thus, the cost in energy in creating an electron-hole 
pair from the Dirac sea is comparable to the Coulomb 
scattering matrix element between pairs of opposite chi- 
rality when there is a single superfluid pair in the system 
(N s = 2). This suggests that creating Dirac electron- hole 
pairs plays an important role in the promotion of pairs of 
one chirality from the bottom to the top of the band. As 
more such promotion of half the pairs from the bottom 
to the top of the band occurs, the wavevector go gradu- 
ally increases and the Coulomb matrix elements decrease 
until it becomes too expensive to create Dirac sea pairs 
and the process stops. The resulting qo determines the 
London penetration depth. 

What happens if the wavefunction of the superconduc- 
tor has such an admixture of states with electron-hole 
pairs created from the Dirac sea, even if very small? The 
new electrons that pop out of the Dirac sea add to the 
electrons in the metal giving rise to a negative charge 
density slightly larger than the compensating positive 
ionic charge density. What does a metal do with excess 
negative charge? It pushes it to the surface! 44J. This 
gives rise to the extra negative charge near the surface 
of superconductors predicted by the theory of hole su- 
perconductivity, shown schematically in Figure 11. The 
compensating positive holes created in the Dirac sea in- 
stead are uniformly distributed over the volume of the 
superconductor, giving rise to the interior positive charge 
density shown in Figure 11. 

For superconductors of linear dimension smaller than 
the London penetration depth charge expulsion does not 
occur [l(J, suggesting that electron- hole pairs in the Dirac 



+ 



.+ 



+ \\ 



^ +\ + / J 



FIG. 11: Schematic picture of a superconducting body (from 
Ref. [lOj]). Negative charge is expelled from the bulk to a sur- 
face layer of thickness Al, the London penetration depth, and 
has charge density p_ given by Eq. (2). As a consequence, 
positive charge density and an outward pointing electric field 
E exist in the interior. 



sea are not being created. This is consistent with the fact 
that the energy scale of Coulomb repulsion at distances 
smaller than 2Al (wavevectors larger than qo) is not large 
enough to create electron-hole pairs in the Dirac sea. For 
such a case the spacing between k-points in the Brillouin 
zone is larger than qo, hence the 'small' Fermi surface in 
Fig. 3 doesn't really exist, and magnetic fields cannot be 
screened. 

It is interesting to note that if the wavefunction of the 
superconductor has terms with electron-hole pairs cre- 
ated from the Dirac sea, the number of electrons is not a 
good quantum number. The form of the BCS wave func- 
tion, that does not have a definite number of electrons, 
should have alerted us to this possibility long ago. 

If this is correct, the key to superconductivity and the 
Meissner effect is the long-range nature of the Coulomb 
repulsion. This was in fact suspected to be the case in 
the early days of superconductivity researchjilj. Then, 
Bohm and Pines [HI allegedly succeeded in removing the 
"troublesome long-range Coulomb interaction" [43| from 
consideration in metal physics. Perhaps they shouldn't 
have. 

Note that the dispersion relation in the center of the 
Brillouin zone becomes Dirac-like (Fig. 8), thus mimick- 
ing the behavior of electrons and holes in Dirac theory. 
The promotion of electrons from k ~ to the top of the 
band mimics the promotion of electrons from negative 
energy states to positive energy states in Dirac theory 
proposed here. 

The energy scale of an electronic energy band is a few 
eV, which is of order of the quantum confinement energy 
of an electron in a distance ao = h 2 /m e e 2 , the Bohr 
radius: 



E, 



band 



2m e aQ 



(44) 
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The Dirac energy scale is the quantum confinement en- 
ergy of an electron in the quantum electron radius Eq. 
(10) 



Dirac 



2m e r 2 q 



2m e c 



(45) 



The Bohr radius and the quantum electron radius are 
related by 



a 
2 



(46) 



with a = e 2 /he the fine structure constant, and the en- 
ergies are related by 

Eband — { — ) 2 E Di rac = { — ) 2 E Dirac (47) 
Tq 2 

The London penetration depth is related to the Bohr 
radius by 



a = - x (2X L ) 



(48) 



when the superfluid density is n s = l/(47rag)[l3], and 
hence the London penetration depth is related to the 
quantum electron radius by 

r q = -xa Q = -x-x (2A L ) (49) 

Thus, the ratio of the band energy scale to the Dirac 
energy scale is 



Eband — {r q qo)EDi 



(50) 



The orbits of electrons in superconductors have radius 
2Al and the electrons in the ground state spin current 
orbit at speed v® (Eq. 4) so their angular momentum 
is 



L 



m e v. 



S(2A L ) 



h 

2 



(51) 



the same as an electron spinning at the speed of light in 
a circle of radius r q . The kinetic energy of a spinning 
electron at speed c can be understood as the quantum 
zero point energy of a particle confined to the quantum 
electron radius r q and yields its rest mass Eq. (45), simi- 
larly the kinetic energy of an electron in the ground state 
spin current in the superconductor can be understood as 
arising from the quantum confinement of a particle in a 
region of linear dimension 2A^ = g^flil] 



4m e 



(52) 



and we have argued in Ref . [40| that this is the condensa- 
tion energy of the superconductor. 

Thus there is an extraordinary parallel between elec- 
trons in superconductors and real electrons: the super- 
conductor is a giant microscope dilating length scales by a 



factor l/(r q q ) = (2/a) 2 , and the electron in the Cooper 
pair of the superconductor orbiting at radius 2Al with 
speed v® is a giant amplified image of the tiny spinning 
electron itself. In a sense (Eq. (49)), the superconduc- 
tor is a 'giant atom'[45( and the atom is a 'giant elec- 
tron', with 'giant' = 2/a. It is natural to expect that 
pair pro duction will play an important role in a giant 
atom[46L |47j. The excess negative charge density near 
the surface p_ (Eq. (2)) is the superfluid charge density 
en s 'diluted' by the expansion factor (r q qo) . 



IX. DISCUSSION 

In the presence of an electric field, electronic energy 
bands will spin-s plit due to the spin-orbit interaction 
(Rashba splitting [2 2]). Since within the theory of hole 
superconductivity a macroscopic electric field exists in 
the interior of superconductors due to negative charge 
expulsion, spin splitting of the bands is a necessary conse- 
quence. However, we have argued that the spin splitting 
that occurs in the superconductor is enormously larger 
(~ 10 5 times) than that implied by the Dirac spin-orbit 
coupling with the macroscopic electric field and is in fact 
driven by the Coulomb interaction between electrons, so 
it corresponds to a 'dynamic generation of spin-orbit cou- 
pling' as proposed by Wu and Zhang 17] . 

Wc have discussed in this paper a new scenario for su- 
perconductivity within the theory of hole superconduc- 
tivity. This scenario unifies many different elements of 
the theory, which we believe makes it compelling. We 
review these elements below. 

(i) The theory predicts that in the transition to su- 
perconductivity electronic orbits expand from a micro- 
scopic radius to a mesoscopic radius 2Al, and that a 
spin current of characteristic wavevector qo = 1/2 Al 
develops [ill l40l]. In the present scenario this is described 
by the small spin-split Fermi surface in Fig. 3 as well as 
the associated spin-split Fermi surfaces near the corners 
of the Brillouin zone. The process by which holes migrate 
from the corners of the Brillouin zone to the bottom of 
the Brillouin zone in the transition to superconductivity 
is consistent with the real-space picture where the orbits 
continuously expandjlll- Note that all electrons in the 
Fermi sea are affected in both processes. 

(ii) In the previous paper in this series it was ar- 
gued that occupation of higher fc-states in the Brillouin 
zone should be associated with negative charge expulsion 
in real space. In the scenario discussed in this paper, this 
occupation of higher /c-states near the corners of the Bril- 
louin zone is directly related to the development of spin 
splitting and a spin current, i.e. the expulsion of electrons 
of one chirality from the center of the Brillouin zone to- 
wards the corners. In Ref. @ we derived electrodynamic 
equations in the charge and spin sectors which showed an 
intimate relationship between charge expulsion and spin 
current development, however the reason for this connec- 
tion was not intuitively clear from the electrodynamics. 
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(iii) The Coulomb matrix element J was shown to drive 
the 'normal' spin-split state in ref . [l6j . thus it is natu- 
ral that it should also play the role in the development 
of the spin current that we had proposed exists in all 
superconductors [ll|, UH. Here we have shown that in- 
deed it plays a key role under the assumptions that spin 
splitting occurs at the center of the Brillouin zone and 
that the normal state Fermi surface is close to the corners 
of the Brillouin zone. 

(iv) We have argued that experiments conclusively 
show that carriers 'undress' from the electron-ion interac- 
tion in the transition to superconductivity and behave as 
free electrons (26|. This is consistent with a Fermi surface 
developing near the center of the Brillouin zone where 
the states are free-electron-like. 

(v) We have given a series of arguments that the 
Coulomb repulsion between electrons will be reduced in 
the scenario discussed here. This is consistent with su- 
perconductivity being driven by the repulsive Coulomb 
interaction as proposed in the theory of hole supercon- 
ductivity rather than by the attractive electron-phonon 
interaction proposed by BCS. The Coulomb scattering 
matrix element V^(q) between Cooper pairs will quite 
generally be attractive (or least repulsive) for q connect- 
ing the center and the corners of the Brillouin zone. To 
take full advantage of this requires Fermi surface parts 
to reside near the center of the zone and near the zone 
corners, as in the scenario discussed here. Spin splitting 
of states with opposite chirality will reduce the Coulomb 
repulsion for Cooper pair scattering across the Fermi sur- 
face and thus the Coulomb repulsion favors spin splitting 
of the Fermi surfaces. The long range Coulomb repulsion 
between electrons deep in the Fermi sea will be reduced 
by promoting electrons of one chirality from the bottom 
to the top of the band, and electron-hole pairs created 
from the Dirac sea will take advantage of it. 

(vi) The facts that superconductors require a nearly 
filled band, that they carry a spin current, and that they 
expel negative charge from the interior, all being pro- 
posed independently in earlier stages of this theory, are 
seen to be naturally related in the scenario discussed here. 
Even the fact that holes exist deep inside the electronic 
Fermi sea was anticipated by the model used in the early 
stages of the theory[2|, |20j (see Eq. (41)), even though 
that prediction of the model was not recognized at the 
time. 

(vii) In the presence of a charge current, the Fermi sur- 
faces in Fig. 3 will be shifted according to the wavevector 
q, with v s — hq/m e the superfluid velocity. The phase 
space for scattering on the small Fermi surface gets pro- 
gressively smaller and vanishes when the edge of the small 
Fermi surface crosses the point k = 0, i.e. for q — q . At 
that point, which is when the carriers with spin current in 
direction opposite to the charge current come to a stop, 
superconductivity should be destroyed. This is consistent 
with the interpretation of the critical current discussed 
in connection with the Spin Meissner effect [ll[ . 

In summary: we have discussed a scenario wherein 



the topology of the Fermi surface (or more precisely its 
equivalent in the superconducting state) changes quali- 
tatively in going from the normal to the superconduct- 
ing state. We propose that every superconductor has a 
hidden hole 'core', a spin-split Fermi surface enclosing 
hole-like carriers in the fc-space region where the lowest 
electron states reside in the normal state (smooth free- 
electron-like bonding states [26]), of radius qo — 1/2 Al. 
This core of undressed hole carriers with definite chirality 
and free-electron-like behavior 'drives' the system and is 
responsible for the remarkable universal physics of super- 
conductors and in particular the Meissner effect as well 
as the Spin Meissner effect. The detailed form of the 
BCS wavefunction, the connection with Dirac physics, 
and experimental consequences will be discussed in fu- 
ture work. The proposed hole-core of superconductors 
may be experimentally detectable by recently developed 
high resolution Laser ARPES spectroscopy 48j . 



Appendix A: next-nearest neighbor matrix elements 

The hopping amplitude between sites i and j for a band 
arising from local s-orbitals fiijc) is 

Uj = - J d 8 r^(r)[-A-V 2 +^i(r)+tr J (r)]^(r) (Al) 

where Ui(r) is the ionic potential. Assuming ifi is an 
atomic orbital it satisfies (-S 2 /(2m e )V 2 + \]j)tpj = £j<Pj, 
and furthermore assuming ifi and ipj are orthogonal Eq. 
(Al) becomes 



U 3 = - I dV;(r)tt(r)¥>i(r) 



We take as ionic potential 



Ui(r) 



Ze 2 



r-Ri 



(A2) 



(A3) 



with Z = 2 the ionic charge, consistent with having a 
charge-neutral system when the band is full (n e = 2) . 
Furthermore we introduce an ionic wave function |49j] 
Xi(r) which is essentially a 6- function at the position of 
the ion. The hopping amplitude Eq. (A2) is then 



tij = 2 J d3 d 3 rV:(r) ^ (r) __ |Xi(r0| 2 (A4) 

On the other hand the correlated hopping matrix element 
(ii\l/r\ij) between these sites is 



(A% = J ^V^(r)^-(r)__ -|^(r')| 2 (A5) 

Comparing Eqs. (A4) and (A5) we conclude that for 
nearest neighbor matrix elements we will have approxi- 
mately 2(At)jj ~ t^, i.e. the effective nearest neighbor 
hopping Eq. (1), t(n e ) — to — n e At becomes indeed very 
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small as the band becomes full (n e — > 2). Furthermore, 
for next- nearest- neighbor sites the fact that |</?i(r')| 2 is 
more spatially extended than |xi(r')| 2 should play a more 
important role since ipj in Eq. (A4) is further away and 
exponentially decaying, hence Eq. (A5) will become rel- 
atively larger compared to (A4) for further than nearest 
neighbor orbitals. This justifies the assumptions made in 
Sect. 5 that the ratio Ai 2 Ao2 is larger than At /to and 



that to — n e At > and £02 — n e At 2 < 0. 
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